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Abstract 


Composite materials have large number of uncertain parameters associ- 
ated with its manufacturing. It is not possible to control these completely. 
This results in variations in the properties of composite materials. For 
accurate modeling, the material properties are, therefore, considered as 
random variables. The influence of randomness in material properties on 
buckling behavior of graphite/epoxy composite laminate has been analyzed 
in the present study with the help of Monte Carlo Simulation and FEM. 
The analysis has been carried out for composite laminates with different 
edge conditions, aspect ratios(ARs), and ply-orientations to random mate- 
rial properties. The study has also been extended to plates with cut-outs 
using NASTRAN software. It has been found that the normalized buck- 
ling response of the composite plates with random materials is similar for 
various aspect ratios, edge conditions and ply-orientations. Further, for a 
given area of the cut-outs, the buckling load for a laminate with elliptical 
cut-out is lower than circular and rectangular cut-outs. The buckling load 
with circular cut-out is the highest. 
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Chapter 1 


Introduction 

The composite materials are formed by the combination of two or more 
materials on a macroscopic scale. The materials used in composites may 
be called as: reinforcing material and parent or matrix material. Struc- 
tures made of such materials are called as composite structures.Composite 
structures can be fabricated to have better performance in engineering ap- 
plications compared to the conventional metallic materials. . Some of the 
properties that can be improved by forming a composite are: stiffness and 
strength to weight ratios, corrosion resistance, thermal characteristics, fa- 
tigue life and wear resistance etc,. These are the reasons for composites 
finding applications in a variety of systems including aircraft and subma- 
rine structures, space structures, automobiles, sports equipment, medical 
prosthetic devices and electronic circuit boards. 

Presently, aircrafts are being manufacturing with a very high percent- 
age of components made from the composite materials. These results in 
reduction in structural weight consequently the payload can be increased 
for a given engine performance. Besides this, the range can be extended 
and operating efficiency improved. Space vehicles having plan forms like 
plates have been designed and built of high strength temperature resistant 
composite materials. Composite materials used for space vehicles provides 
structural efficiency, ease of fabrication and freedom from thermal distor- 
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tion. 

All materials will have dispersion in material properties due to lack of 
complete control over manufacturing / fabrication / processing techniques 
employed. Composite materials experience larger dispersions in their ma- 
terial properties compared to conventional materials due to larger number 
of parameters involved in their fabrication. Also, the uncertainties involved 
in manufacturing and processing techniques employed for composites are 
more as compared to isotropic materials. These uncertainties can occur 
due to air entrapment, de-lamination, lack of resin, incomplete curing of 
resin, excess resin between layers etc,. Moreover, the different test meth- 
ods provide variable property data. Zweben [1] illustrated the dispersion in 
test data as a result of conventional testing methods used for composites. 

Considering the above aspects, the properties of structural materials in 
general and composites in particular, may be modeled as a random process. 
However, in conventional structural analysis generally used for composites, 
the material properties like elastic modulus, poisson’s ratio etc. are com- 
monly assumed to be deterministic quantities. Studies made by various 
researchers have shown that the effects of material property randomness 
on dynamic response are much more pronounced in case of composites 
as compared to conventional materials [2]. It is, therefore, important to 
consider the material properties as random in the case of composites vis- 
a-vis the metallic materials as otherwise the predicted response may differ 
significantly from the observed values and the structures may not be safe. 

In spite of the well-developed theory and computational methods for 
structural responses for deterministic material properties, the random struc- 
tural analysis where the material properties are considered random remains 
underdeveloped. Nevertheless, for reliability of design and considering the 
sensitivity of the application, accurate predictions of the system behavior 
of structure made up of composites in the presence of the uncertainties 
favors a probabilistic analysis approach for composites by modeling their 
properties as random variables (RVs). 
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1.1 Literature survey 

The uncertainties in strength and stiffness properties of composites has 
been studied by many researchers. Ibrahim [3] reviewed number of topics 
on structural dynamics with parameter uncertainties. Nakagiri et al.[4] 
have studied simply supported (SS) graphite/epoxy plates with Stochastic 
Finite Element Method (SFEM) taking fiber orientation, layer thickness 
and number of layers as random variables, and found that the overall stiff- 
ness of Fiber Reinforced Plastic(FRP) laminated plates is found out to be 
largely dependent on fiber orientation. Leissa [5] and Martin have ana- 
lyzed the vibration and buckling of rectangular composite plates and have 
established that variation in fiber spacing or redistribution of fibers tend 
to increase the buckling load by 38 percent and the fundamental frequency 
by 21 percent. Englestad and Reddy [6] studied metal matrix compos- 
ites (MMC) based on probabilistic micro mechanics nonlinear analysis. 
They have used Monte Carlo Simulation(MCS), and different probabilistic 
distributions to incorporate the uncertainty in basic material properties. 
Vinckenroy and de Wilde [7] in their first part of the work have established 
a procedure to obtain the best fit for each material property. Further, they 
have studied the behavior of perforated plate and determined the probabil- 
ity distribution of the response. The input variables have been simulated 
using MCS and SFEM. Salim et.al [8] studied the statistical response of the 
system considering the material properties such as, longitudinal modulus, 
poisson’s ratio, transverse elastic modulus, etc., as independent random 
variables on [90°/0°/90°] and 0°, rectangular, SS, and CC plates subject to 
transverse static loads. Salim et.al [9] have shown that there is change in 
the rate and the order of the natural frequency with SD of input RV and 
boundary conditions of the plate. Raj et.al [10] studied the response of 
composite plates with random material properties using FEM and MCS, 
and found that the longitudinal modulus En and in-plane shear modu- 
lus (?n are the most critical material properties influencing the deflection 



4 


characteristics of the laminate. Jagdeesh [11] studied the static response 
of composite laminates with randomness in material properties as well as 
external loading. There appears to be a very few published literature 
when both material properties and external loading are modeled as ran- 
dom variables, especially in the area of composites. Salim et.al [12] used 
the perturbation technique with Rayleigh-Ritz (RR) formulation to ana- 
lyze the bending, buckling and vibration of composite plates. The RR 
method can be used only for regular boundary problems. This problem 
can be over come by using FEM for buckling analysis. Though the per- 
turbation technique though gives acceptable results, the method fails for 
larger dispersions in the input random variables. The SFEM can handle 
larger uncertainties in the input parameters. However, it is computation- 
ally expensive and being a specialized technique its other applications are 
limited. 


1.2 Present work 

With the usage of fiber reinforced and laminated composite materials, the 
domain of application and range of various shapes of plan forms like plates 
has enormously increased. Aircraft as well as space vehicles, having plan 
forms like plates, have been designed and successfully built of high strength 
temperature resistant composite materials. The skin of aircraft structures 
are also composed of plate forms, built of stiffened plates and /or compos- 
ite materials. Thus, with such immense potential of applications in various 
fields of engineering, the design criteria becomes stringent. Therefore, it 
is of paramount importance to consider the basic material parameters as 
random and study their effect on the response characteristics of plate struc- 
tures. The difficulties in applications of perturbation technique and SFEM 
to random material properties in problem of composite laminates has moti- 
vated the use of a FEM and commercially available NASTRAN software in 
the present study. The material properties are assumed to vary randomly. 
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The stiffness properties of lamina such as E n ,E X 2 ,G l2 and u i2 have been 
selected to represent its behavior and are taken to be random variables. 
In the present analysis the buckling of composite laminates due to uniaxial 
in-plane load has been studied with the above material properties modeled 
as random variables. MCS and FEM is adopted to obtain the critical load 
of the laminated plates and its statistics with different input characteristics 
and cut-out shapes. 


1.3 Layout of the thesis 

In the present work, the buckling characteristics of the plate with random 
material properties have been studied. The thesis is organized in the fol- 
lowing way. The formulation for buckling analysis has been presented in 
Chapter 2. The Monte Carlo simulation has been used for random anal- 
ysis. This has been discussed in Chapter 3. The response characteristics 
of the plate have been presented in Chapter 4. The conclusions and a few 
suggestions for further work are given in Chapter 5. 



Chapter 2 


Formulation for Buckling 


2.1 Stress-Strain relations for orthotropic lamina 


The constituent equations for kth orthotropic lamina of the laminate in a 
two-dimensional state of stress along 1 and 2 Can be written as 



Q 12 
Q22 

0 


0 

0 

Qg 6 



( 2 . 1 ) 


It can also be written as 



where 


lQij] ilj=1 , 2 ,...,e are the stiffness coefficients of the lamina and they are related 
to the material properties as follows 


fi 



7 





W 


Figure 2.1: Laminated composite rectangular plate 
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Stress-strain relations along x and y axes for the lamina with fiber orien- 
tation 9 (as shown in Fig. 2.1) are 
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where 


are the transformed stiffness coefficients for the kth lamina 


[QiJ 

The relations between the elements of the |Q matrix and the [Q] are [13] 


Qn — Qnn 4 + Q22™ 4 + 2 (Qn + 2 Qe$) m 2 n 2 
Q22 = Qn ^ 4 + Q22E 4 + 2 (Q12 + 2Qqq) m 2 n 2 
Q12 = (Qn + Q22 — 4 Q 6 e) m 2 n 2 + Q i2 ( rn 4 + n 4 ) 

Q 66 = (Q11 + Q22 — 2Q12 — 2 Q 66 ) m 2 n 2 + Qee (m 4 + n 4 ) 
Qi 6 = (Qn - Q12 - 2 Q 66 ) mn 3 - (Q 22 - Q12 - 2 Q 6 e) nm 3 
Q26 = (Qn - Q12 - 2 Q 66 ) nm 3 - (Q 22 - Q12 - 2 Q 6 s) mn 3 
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in which 

m = sin# and n = cos 9. 

Considering only the transverse displacement w(x,y), and symmetric lam- 
inates, the strain displacement relations are given by 

We note, that, for antisymmetric angle-ply laminate, the effect of the 
coupling terms B 13 and B 2 3 are reduced when the 9° and -9° plies are inter- 
spersed like [9°/ — 9°/9 0 / — 9°]. Hence, the number of 9° and —9° plies in an 
angle-ply antisymmetric laminate are increasingly interspersed, the lami- 
nate behaves as if were predominantly symmetric, and having membrane 
and bending orthotropy, since the effect of the coupling terms (.4 i 6 =0, 
.4 26 =0) and no bend twist terms (D I6 = 0, D 26 =0) present for any any anti- 
symmetric laminate configuration. 



Therefore 

the stress-strain relations for the kth lamina are 
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It can also be written as 
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called as bending or flexure stiffness matrix. 


2.3 Finite Element Formulation 

For a two-dimensional state of stress the strain energy stored in the kth 
lamina can be expressed as 



From Eq. (2.7) 

G x ® xy Txy 

Substituting Eq.(2.12) in Eq.(2.11) we can get 
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where {D 2 w} = < 




_ d 2 ^ 

dx 2 

0 2 w 
dy 2 

dxdy 


(2.15) 


Integrating over the whole thickness of the laminate, we obtain the total 
strain energy stored 


U = j A {^y [FS]{D 2 w}clA (2.1G) 

where 

[i ? 5]= bending stiffness matrix =[D;j]. Work done by the external force 



Therefore the variational functional for the present problem 


i=u ~ w= \L [{■ £)2 ”'} r |i?51 


dA 


(2.17) 


(2.18) 


Now 


w = {N} eT {5} e 

= {5} eT {TV 6 } (2.19) 


where 

(iV} e is the shape function vector and {5} e is vector of elemental degrees of 
freedom. 
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The shape functions used for a four noded rectangular plate element with 
four nodal d.o.f (w, 9 X , 9 y , 9 xy ) are 

Ni = h (0 h (17) 

N 2 = a 7i (0 fi (rj) 

W = * 7 i (0 9i (rj) 

^ = a e b e gx (fl ^ (77) 

= / 2 (0 fi (V) 
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= *72 (0 * ( 77 ) 

iVg = a e b e g2 (£) pi (77) 

= h (0 h ( V ) 

iVfo = °72 (0 / 2 (77) 

= 672 (0 P2 (77) 

7 Vf 2 = a e i e p 2 (0 p 2 (77) 

^3 = /1 (0 h (7?) 

Ni £ 4 = a e p! (0 / 2 (77) 

^ = Vfi (0 92 (77) 

^6 = « e 6 e p! (0 92 (77) 

/1 (s) = J (2 - 3s + s 3 ) 

/ 2 (s) = — ^2 + 3s + s 3 j 

( s ) = ^ (l — s _ S 2 + ^ 3 ) 

52 ( 5 ) = ^'(-1 - S + s 2 + s 3 ) 
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{r\= ( N x {B x } e {B x Y T dA (2.24) 

J A c 

where [k c ] is the elemental stiffness matrix and [ f e ] is the elemental force 
matrix. Substituting Eq.2.21 and Eq.2.22 in Eq.(2.18) we can get 

I=\J A [{<S}‘ T [cf [re] [Cl* {{}' - N, {if {B,Y {B,f {«}'] (2.25) 

Substituting Eqs.(2.23) ahd (2.24) in the above equation, we get 

/ = 5 £ [{n' r [*]‘ W* - {if [/]• {*}'] (2-26) 

1 e=l 

The global stiffness matrix, global force matrix and vector global degrees 
of freedom can be obtained as 


[ k \ = EM' P- 27 ) 

e—l 

ne 

[F] = £ [/)• (2.28) 

e=l 

ne 

[W] = £ {5}* (2.29) 

e—l 

where [ K] e are obtained from [k] e using the following simplified assembly 
equations. 

^4i — 3,4s-3 = ^4p-3,4?-3 
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Otherwise [iiT] e =0 

Similar relations are valid for obtaining the [F] e from [f] e for determining 
{H^} e , the following relations can be used. 
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when r = C ep . 

Otherwise {iy} e = 0 

where C ep is the element of the connectivity 

r , I" 1 2 3 4 

l°}= .... ( 2 - 31 ) 

• • • • 

• * • ♦ 

where 

the number of rows= number of elements and number of columns =number 
of nodes per element. In (f , 77 ) coordinates 

[C'l = [t] [C'l (2.32) 

where 

(I) 2 0 0 

0 (|) 2 0 (2.33) 

0 ° Sg. 

This has been obtained by using the chain rule of differentiation and using 
the given conversion for rectangular elements from x, y to natural coordi- 

9x Q 

nates £, rj Jacobian matrix [J] = 

0 & 

L 

Det [J] = a e .b e 
x = a 1 + a e f; 
y = h + b e r] 

Substituting Eqs.(2.27),(2.28) and (2.29) in Eq.(2.26) we get 

/ = i pvf \K\ {W} - {Wf [F] {»'} (2.34) 

Minimizing the above variational functional 



d{W} 0 


(2.35) 
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From the above equation, we get 

[K]{W}-[F]{\V} = 0 (2.36) 

Thus the final FEM equation for the present problem can be written as 

WOn-m{^} = 0 (2.37) 

This is an eigen value problem with characteristic equation \[K] - [F]| = 
0, which has been solved to obtain the buckling loads and the mode shapes 
by using NAG routine. 

2.4 Buckling analysis using NASTRAN 

In the present problem the major analysis is carried out by NASTRAN 
finite element method analysis. At buckling the determinant of the sum of 
the elastic stiffness matrix and differential stiffness matrix is equal to zero. 

K e + K%\ = 0 (2.38) 

Where K e is the elastic stiffness and K is the differential stiffness matrix 
at critical buckling load 

Ki (u) = A x K d (u) (2.39) 

K e + XK d {u)\ = 0 (2.40) 

The above equations is known as the eigen value problem with eigen values 
Ai, A 2 , • • • , A n where N is the size of the K e and K d matrices (and is equal 
to the degrees of freedom in the analysis set of the problem). Each value 
of A corresponds to a buckling load. 

2.5 Methodology for cut-outs 

Cut-outs are present in laminated plates in many applications. In the 
present study three different shapes of cut-outs have been taken namely, 
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rectangular, circular, and elliptical. For comparison of the plate behavior 
all the cut-outs have been taken to have the same area. 

1. For rectangular cut-out as shown in Fig. 3.1 
A c = c x d 

Where 

A c =area of cut-out 
c=length along major axis 
d=length along minor axis 

2. For circular cut-out as shown in Fig. 3.1 

A c = 7rr 2 

r=radius of the cut-out 

3. For elliptical cut-out as shown in Fig. 3.1 

A c = ircdf4- 



Chapter 3 


Solution approach-Monte Carlo 
simulation technique 


The problem of analysis of composite structures with variation in material 
properties using Monte Carlo simulation approach can be divided in to the 
following steps. 

1. Random input generation. 

2. System analysis. 

3. Post processing of response behavior. 


3.1 Random input generation 

3.1.1 Number generation 

A set of random numbers are generated of a given sample size, mean, stan- 
dard deviation and probability distribution. Standard library subroutines 
are available for this. These can be known as pseudo random numbers as 
these are artificially generated through numerical subroutines. 
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3.1.2 Processing of raw data 

For getting the random numbers with desired characteristics, the pseudo 
r an dom numbers generated by G05FDF from NAG library are prepro- 
cessed before feeding them as input to the FEM analysis. The samples 
are checked for shift in the input mean and standard deviation which are 
usually present in a small amount because of finite length arithmetic of 
computers. For adjusting these shift these three Arithmetic formulations 
can be used. 

1. Modification of the random numbers to shift the mean. 

TlTi =Xi- [fi e - fir ) (3-1) 


2. Modification of the random number to change the standard deviation 
(SD) 

nr . 

(3.2) 


Xi 

nr* = — x o T 

Oe 


3. Modification of the random number to change the SD and shift the 
mean simultaneously. 

[Xi - fi e ) X <7 r 


nTi — 


+ fir 


(3-3) 


where n 

X>i 

fie = 

n 

E - v ) 2 

<?e = \ — T 

\ n — 1 

H r = required mean values 

cr r = required standard deviation values 

nrj=new random numbers. 


(3.4) 


(3.5) 
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3.2 System analysis 

The governing equation of the buckling analysis of laminated composite 
has been solved with the help of finite element technique. A 4-DOF C l 
continuity model with 4-noded Lagrangian element has been employed in 
the present problem. In the finite element program the various stages are. 

1. pre-processor 

2. processor, and 

3. post-processor 

3.2.1 pre-processor 

It supplies the following data required by the processor 

1. Mesh data 

• number of elements 

• number of nodes per element 

• number of DOF per node 

• connectivity matrix 

• vector of nodal co-ordinates 

2. vector of nodal values of geometric properties, material properties and 
other properties. 

3. specifying the node numbers and the designated boundary values. 

3.2.2 processor 

Finite element equations are assembled and solved. 

• calculation of elemental coefficient matrix and right side vector 


• global assembly 
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• application of essential boundary conditions 

• solving the equations 

3.2.3 post-processor 

This involves the calculating and plotting of required results. The sample 
of the buckling response of the plate is obtained from the FEM analyzer 
and NASTRAN software. Its statistical information is extracted from the 
Eqns. (3.4) and (3.5). 






Chapter 4 


Results and Discussions 

4.1 Convergence study 

Convergence study has been carried out for three types of plates namely 
isotropic square plate, composite square plate and composite square plate 
with circular cut-out. Table 4.1(A) shows the convergence for square 
isotropic plate with material properties E — 200 GPa and u = 0.28. The 
convergence study for composite square plate and composite square plate 
with circular cut-out as shown in Table 4.1(B) and (C) for material prop- 
erties E n = 206.8425GPa, E n = 5.171 GPa, u n = 0.25 and G l2 = 2.585 GPa. 
Typical results for a circular cut-out with a/d — 5.0 are shown in Table 
4.1(C). It it observed that a 6x6 mesh gives sufficient convergence from en- 
gineering point. Hence, further computations have been carried out based 
on 6x6 mesh for the full plate. 

4.2 Validation study 

Validation studies have been carried for symmetric square cross ply-laminate 
with simply supported boundary conditions for different mesh sizes. These 
results are presented in Table 4.2. The material properties used for this 
are also shown in Table 4.2 

From the Table 4.3 it is observed that for all edges simply supported 
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and all edges clamped for a given minor axis, as the major axis reduces the 
buckling load increases and then reduces at certain value of a/c which is 
close to given a/d ratio. Similar observations can be made for a plate which 
is clamped on two edges and two edges simply supported. In the case of a 
plate with three edges clamped and one side free, the trend is reverses. 
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Table 4.1: Convergence study for square isotropic and composite plates 


EDGE CONDITIONS 

MESH SIZE 

4 

CRITICAL LOAD! 10 >N 


C 


2X2 

158.5388 

C 


c 

4X4 

159.4702 

6X6 

159.5098 


C 



s 


2X2 

63.2396 

S 


s 

4X4 

64.2293 




6X6 

64.2691 


s 





( A ) : Convergency study for Isotropic square plate 


EDGE CONDITIONS 

MESH SIZE 

4 

CRITICAL LOAD (10 ) N 


c 


2X2 

79.13348 

C 


c 

4X4 

30.58752 

6X6 



30.5831 5 


c 



s 


2X2 

7.932761 

S 


S 

4X4 

12.40192 




6X6 

12.40335 


S 



( B ) : Convergency study for Composite square plate 


EDGE CONDITIONS 

MESH SIZE 

4 

CRITICAL LOAD ( 10 ) N 



c 


2X2 

19.705526 

1 

r 



£ 

c 

o 

[<* d *4 

c 

4X4 

21.386021 

u 


6X6 

21.391097 



c 




s 


2X2 

5.2786344 






s 


s 

4X4 

8.0285583 





6X6 

8.0310422 



s 



( C ) : Convergecy study square composite plate with circular cut-out (a/d=5) 
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Table 4.2: Validation of critical loads for a square simply supported composite plate 


Laminate(deg) 

Reference 

Critical loads ( N ) 

Number of eiemets in the square plate 

[14] 

2x2 

4x4 

6x6 

8x8 

[0/901. 

18.20341 

18.2232534 

18.2416366 

18.2076981 

18.206495 

[0/90/0] 

s 

18.19803 

18.1270658 

18.20268707 

18.1991125 

18.202405 

[0/90/0/90] 

s 

18.19803 

18.0794632 

18.20167619 

18.1987803 

18.196925 

[0/90/0/90/0^ 

18.20098 

18.0387587 

18.20086620 

18.1977916 

18.197267 


Ejj= 140 GPa, E |2 =10 GPa, Gj 2 = S GPa and V |2=0J 


Table 4.3: Critical load for elliptic cut-outs with different boundary conditions for 
deterministic symmetric composite laminate 


Critical loads 

geometry 

edge conditions 
Dotations 1X34 

a /c=4 

a/c=5 

a/c=6 

a/c=7 

a/c=8 




cccc 

5X70304 

5X9248 

5.6126 

5.74S87 

5381212 

1 

c± 

2 

a 

2.66907 

2.65275 

2.55812 

236623 

238523 

CCCF 


4 

a 







a/d=10 


CCFF 

1.982127 

2.05487 

1.942303 

1.88537 

1.97432 




CSCS 

3320271 

3.4449 

3.70915 

3.8545 

3.79850 




SSSS 

2.038672 

2.08123 

2.1298 

2.2038 

2.09114 
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4.3 Plates with random material properties 

The approach presented in the previous chapters is used to find the buckling 
load statistics for rectangular composite plates. In section 4.4 behavior of 
rectangular plate with different boundary conditions and ply-orientations 
has been presented assuming the material properties as random. In sections 
4.5 - 4.7 results have been presented for the plates with rectangular, circular 
and elliptical cut-outs respectively. Initially the effect of randomness in 
individual material property has been investigated to identify the most 
critical parameter. Subsequently only two parameters En and G n have 
been taken as random variables. 

• Samples are generated with the standard deviation varying from 5 
percent to a 20 percent of their mean values with sample-size 1000 
and 10 such samples are used for each study. 

• The material properties considered to be random are E n , E x2 ,vn and 
G 12 * 

Results have been computed for graphite/epoxy composites with mean 
values of the material properties as as given below [14]. 

E n = 206.8425G.Pa 
E n = 5.17lGPa 
un = 0.25 and 
G n = 2.585 GPa 


4.4 Buckling characteristics of square and rectan- 
gular plates 

In this section a study of the effect of dispersion in single input random 
variable while the characteristics of the other input variables are kept con- 
stant on plates with different aspect-ratios (AR), boundary conditions and 
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ply-orientations is carried out. Two classical boundary conditions have 
been considered namely all edges SS and all edges CC. The critical loads 
are normalized with respect to the critical loads for[0°/90 0 /0°/90 0 ] lami- 
nate treating the material properties to be deterministic. Three different 
lay-ups namely [45°/— 45°/— 45°/45°] , [45°/— 45°/45 0 /— 45°] and [0°/90°/0° /90°] 
laminates have been studied. 

Figs. 4.1 to Figs. 4.8 show the mean and standard deviation of normal- 
ized critical load corresponding to each input variable as random . Each 
figure represents six different studies. These presents the effect of the SD 
of input RV on the buckling load of the plates with different boundary 
conditions, ARs and ply-orientations. 

The Figs. 4.9 and 4.10 show the mean and standard deviation of nor- 
malized critical load when both E n and G n are random variables. 

4.4.1 Effect of longitudinal elastic modulus En 

Fig 4.1 shows the variation of mean of the buckling load with standard 
deviation of E n . While Fig 4.2 shows the variation of standard deviation 
of buckling load with standard deviation of E n . The response shows a 
generally a linear behavior. The plate with all edges CC is affected more 
by a change in input randomness compared to a plate with all edges SS. 
The plate shows least sensitiveness to input dispersion as the aspect ratio 
increases. From the limited study conducted it is observed that cross-ply 
laminate is less sensitive compared to angle-ply laminate for all edges SS 
due to change in a / n of En . 

4.4.2 Effect of transverse elastic modulus Eu 

Fig 4.3 shows the variation of mean of the buckling load with standard 
deviation of E 12 • While Fig 4.4 shows the variation of standard deviation 
of buckling load with standard deviation of En- The general behavior is 
almost linear, but unlike the case £ u > here the slope of both mean and 
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standard deviation of normalized critical load curves are decreasing trend. 
The effect of random variable E n is less compared to random variable E n . 
As observed in the case of En here also the plate is less sensitive to input 
dispersion as aspect ratio increases. All edges CC plate is more affected by 
dispersion in E\ 2 than all edges SS. The cross-ply laminate is least affected 
compared to angle-ply laminate for all edges SS due to change in a/n of 
En- 

4.4.3 Effect of major Poisson’s ratio v \2 

Fig 4.5 shows the variation of mean of the buckling load with standard 
deviation of vn- While Fig 4.6 shows the variation of standard deviation 
of buckling load with standard deviation of v i2 - 

The general nature of the curves are mostly linear. The plate with all 
edges CC is more affected than plate with all edges SS. From the com- 
parison of all variables it has been found that the random variable i/ 12 has 
less effect compared to random variables E n , En and Gn- Here too the 
cross-ply laminate is least affected than angle-ply laminate for all edges SS 
due to change in aj\x of vn- 

4.4.4 Effect of rigidity modulus G 12 

Fig 4.7 shows the variation of mean of the buckling load with standard 
deviation of Gn- While Fig 4.8 shows the' variation of standard deviation 
of buckling load with standard deviation of Gn - 

Here too the characteristics are mostly linear. The plate with all edges 
CC is more affected than plate with all edges SS. Here the slope of re- 
sponse characteristic curves have an increasing trend compared to En and 
Vn response characteristic curves. Next to the random variable En, the 
random variable Gn has been found to affect the buckling characteristics. 
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4.4.5 Combined effect of longitudinal modulus Eu and rigid- 
ity modulus G 12 

Fig 4.9 shows the variation of mean of the buckling load with standard 
deviation of Eu and G n- While Fig 4.10 shows the variation of standard 
deviation of buckling load with standard deviation of E n and G 12 * 

The [0 o /90 o /O°/9O°] laminate is least affected compared to [45°/ — 45°/ — 45°/45°] 
and [45°/— 45°/45 0 / — 45°] laminates with all edges SS. The plate with all 
edges CC is more affected than plate with all edges SS. The general be- 
havior is almost linear, unlike the cases of random variables E Ul Eu, vm 
Gui here the slope of mean and standard deviation curves are increasing 
trend. 
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Figure 4.4: Plate buckling(SD)characteristics with En as random 
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Figure 4.5: Plate buckling(Mean) characteristics with u l2 as random 
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4.5 Buckling characteristics of plates with rectan- 
gular cut-out. 

In aerospace structure applications weight reduction is an important design 
criteria. Besides being required due to functional requirements, cut-out 
also have relevance to the problem of weight reduction. More and more 
structures are used with cut-outs. This motivated the study of the behav- 
ior of composite plates with cut-outs. A standard plate has been taken 
with AR=1. Here the cut-out has been used in rectangular plates with 
cut-out side ratio (c/d) having values 1.15, 0.96 and 0.87. These values 
have been chosen such that the areas for different shapes of cut-out is con- 
stant. The normalization for the mean and SD of critical loads have been 
done with respect to a cross-ply laminate without cut-out. Results pre- 
sented graphically to bring out the effect of boundary conditions, ARs and 
ply-orientation for various c/d ratios mentioned earlier. Three different 
lay-ups namely [45°/-45°/-45 0 /45 0 ], [45°/ -45°/45°/-45°] and [0 o /90°/0°/90 o ] 
laminates have been studied. 

Figs. 4.11 to Figs. 4.18 show the mean and standard deviation of 
normalized critical load corresponding to each input variable as random . 
Each figure represents six different studies. These present the effect of the 
SD of input RV on the buckling load of the plates with different boundary 
conditions, ARs and ply-orientations. 

The Figs. 4.19 and 4.20 show the mean and standard deviation of nor- 
malized critical load when both E\\ and G 12 are random variables. 


4.5.1 Effect of longitudinal elastic modulus En 

Fig 4.11 shows the variation of mean of the buckling load with standard 
deviation of En. Fig 4.12 shows the variation of standard deviation of buck- 
ling load with standard deviation of En. The response shows a generally 
linear behavior. The plate with all edges CC is more affected by a change in 
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input randomness compared to a plate with all edges SS. The plate shows 
least sensitiveness to input dispersion as the aspect ratio increases. From 
an limited study conducted it is observed that cross-ply laminate is less 
sensitive compared to angle-ply laminate for all edges SS due to change in 
ajp of Eu. 


4.5.2 Effect of transverse elastic modulus Eu 

Fig 4.13 shows the variation of mean of the buckling load with standard 
deviation of E u . Fig 4.14 shows the variation of standard deviation of 
buckling load with standard deviation of En- The general behavior is 
almost linear, but unlike the case En, here the slope of both mean and 
standard deviation of normalized critical load curves are decreasing trend. 
The effect of random variable E i2 is less compared to random variable Eu. 
As observed in the case of Eu here also the plate is less sensitive to input 
dispersion as aspect ratio increases. All edges CC plate is more affected by 
dispersion in E i2 than all edges SS. The cross-ply laminate is least affected 
compared to angle-ply laminate for all edges SS due to change in cr/p of 
E 12 - 

4.5.3 Effect of major Poisson’s ratio uu 

Fig 4.15 shows the variation of mean of the buckling load with standard 
deviation of u i2 . Fig 4.16 shows the variation of standard deviation of 
buckling load with standard deviation of iy i2 . 

The general nature of the purves are mostly linear. The plate with all 
edges CC is more affected than plate with all edges SS. From the com- 
parison of all variables it has been found that the random variable u n has 
less effect compared to random variables En, E\ 2 and Gi 2 . Here too the 
cross-ply laminate is least affected than angle-ply laminate for all edges SS 
due to change in cr//x of u v2 . 
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4.5.4 Effect of rigidity modulus G u 

Fig 4.17 shows the variation of mean of the buckling load with standard 
deviation of G 12 - Fig 4.18 shows the variation of standard deviation of 
buckling load with standard deviation of G 12 . 

Here too the characteristics are mostly linear. The plate with all edges 
CC is more affected than plate with all edges SS. Here the slope of re- 
sponse characteristic curves have an increasing trend compared to E u and 
uu response characteristic curves. Next to the random variable E n , the 
random variable G n has been affected more for buckling characteristics. 

4.5.5 Combined effect of longitudinal modulus E n and rigid- 
ity modulus G 12 

Fig 4.19 shows the variation of mean of the buckling load with standard 
deviation of E n and Gi 2 . Fig 4.20 shows the variation of standard deviation 
of buckling load with standard deviation of En and G 12 . 

The [0°/ 90 o /0°/90°] laminate is least affected compared to [45°/ —45°/ — 45°/45°] 
and [45°/ — 45°/ 45°/ — 45°] laminates with all edges SS. The plate with all 
edges CC is more affected than plate with all edges SS. The general be- 
havior is almost linear, unlike the cases of random variables En, En, v X2 , 

G 12 , here the slope of mean and standard deviation curves are increasing 
trend. 
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Figure 4.13: Buckling(Mean) characteristics of plate with rectangular cut-out for E \ 2 
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Figure 4.17: Buckling(Mean) characteristics of plate with rectangular cut-out for G 
as random 
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Figure 4.18: Buckling(SD) characteristics of plate with rectangular cut-out for G\2 as 


random 
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4.6 Buckling characteristics of plates with circular 
cut-out 

In this section the circular cutouts has been used in rectangular plates with 
plate length to cut-out ratio a/d having 1.58. A standard plate is taken with 
AR = 1. The normalization for the mean and SD of critical loads have 
been done with respect to a cross-ply laminate without cut-out. Results 
presented graphically to bring out the effect of boundary conditions, ARs 
and ply-orientation for various a/d ratio mentioned earlier. Three diff erent 
lay-ups namely [45°/— 45°/— 45°/45°], [45°/— 45°/45 0 /— 45°] and [0°/90°/0 o /90 o ] 
laminates have been studied. 

The Figs. 4.21 to Figs. 4.28 show the mean and standard deviation of 
normalized critical load corresponding to each input variable as random . 
Each figure represents six different studies. These presents the effect of the 
SD of input RV on the buckling load of the plates with different boundary 
conditions, ARs and ply-orientations. 

The Figs. 4.29 and 4.30 show the mean and standard deviation of nor- 
malized critical load when both En and G 12 are random variables. 

4.6.1 Effect of longitudinal elastic modulus E\\ 

Fig 4.21 shows the variation of mean of the buckling load with standard 
deviation of E\\. Fig 4.22 shows the variation of standard deviation of buck- 
ling load with standard deviation of E n . The response shows a generally 
linear behavior. The plate with all edges CC is more affected by a change in 
input randomness compared to a plate with all edges SS. The plate shows 
least sensitiveness to input dispersion as the aspect ratio increases. From 
an limited study conducted it is observed that cross-ply laminate is less 
sensitive compared to angle-ply laminate for all edges SS due to change in 
a jn of En. 
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4.6.2 Effect of transverse elastic modulus E u 

Fig 4.23 shows the variation of mean of the buckling load with standard 
deviation of E\ 2 . Fig 4.24 shows the variation of standard deviation of 
buckling load with standard deviation of E n . The general behavior is 
almost linear, but unlike the case E n , here the slope of both mean and 
standard deviation of normalized critical load curves are decreasing trend. 
The effect of random variable E 12 is less compared to random variable 
En . As observed in the case of En here also the plate is less sensitive to 
input dispersion as aspect ratio increases. All edges clamped plate is more 
affected by dispersion in E l2 than all edges SS. The cross-ply laminate is 
least affected compared to angle-ply laminate for all edges SS due to change 
in a/fi of E 12 . 

4.6.3 Effect of major Poisson’s ratio vyi 

Fig 4.25 shows the variation of mean of the buckling load with standard 
deviation of u l2 . Fig 4.26 shows the variation of standard deviation of 
buckling load with standard deviation of i/ 12 . 

The general nature of the curves are mostly linear. The plate with all 
edges CC is more affected than plate with all edges SS. From the com- 
parison of all variables it has been found that the random variable i/ 12 has 
lower effect compared to random variables E n , E i2 and G 12 . Here too the 
cross-ply laminate is least affected than angle-ply laminate for all edges SS 
due to change in (j/jx of i/ 12 . 

4.6.4 Effect of rigidity modulus G 12 

Fig 4.27 shows the variation of mean of the buckling load with standard 
deviation of G 12 . Fig 4.28 shows the variation of standard deviation of 
buckling load with standard deviation of G \ 2 . 

Here too the characteristics are mostly linear. The plate with all edges 
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sponse characteristic curves have an increasing trend compared to E n and 
U 12 response characteristic curves. Next to the random variable E n , the 
random variable G n has been more affected for buckling characteristics. 

4.6.5 Combined effect of longitudinal modulus E n and rigid- 
ity modulus G 12 

Fig 4.29 shows the variation of mean of the buckling load with standard 
deviation of En and G 12 . Fig 4.30 shows the variation of standard deviation 
of buckling load with standard deviation of E u and (? 12 . 

The [0° / 90° / 0° / 90°] laminate is least affected compared to [45°/— 45°/— 45°/45°] 
and [45°/— 45°/ 45°/— 45°] laminates with all edges SS. The plate with all 
edges CC is more affected than plate with all edges SS. The general be- 
havior is almost linear, unlike the cases of random variables E n , E l2 , u x2 , 

G\ 2 i here the slope of mean and standard deviation curves are increasing 
trend. 
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Figure 4.21: Buckling(Mean) characteristics of plate with circular cut-out for En as 


random 
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Figure 4.22: Buckling(SD) characteristics of plate with circular cut-out for E n as 


random 
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Figure 4 . 25 : Buckling(SD) characteristics of plate with circular cut-out for ^12 as 


random 
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Figure 4.27: Buckling(Mean) characteristics of plate with circular cut-out for G 12 as 


random 
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Figure 4.30: Buckling(SD) characteristics of plate with circular cut-out for En and 
G 12 as random 
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4.7 Buckling characteristics of plates with elliptical 
cut-out 

Results have been obtained for elliptical cut-out in rectangular plates with 
cut-out ratios (c/d) having values 2.5, 1.225 and 0.4. The area of cut-out 
has been kept same as that of the circular and rectangular cut-outs. A 
standard plate is taken with AR = 1. The normalization for the mean 
and SD of critical loads have been done with respect to a cross-ply lam- 
inate without cut-out. Results presented graphically to bring out the ef- 
fect of boundary conditions, ARs and ply-orientation for various c/d ra- 
tios mentioned earlier. Three different lay-ups namely [45°/— 45°/— 45°/45°], 
[45°/— 45°/45 0 /— 45°] and [0°/90 o /0°/90 o ] laminates have been studied. 

The Figs. 4.31 to Figs. 4.38 show the mean and standard deviation of 
normalized critical load corresponding to each input variable as random . 
Each figure represents six different studies. These presents the effect of the 
SD of input RV on the buckling load of the plates with different boundary 
conditions, ARs and ply-orientations. 

The Figs. 4.39 and 4.40 show the mean and standard deviation of nor- 
malized critical load when both E n and Gn are random variables. 

4.7.1 Effect of longitudinal elastic modulus E\\ 

Fig 4.31 shows the variation of mean of the buckling load with standard 
deviation of Eu. Fig 4.32 shows the variation of standard deviation of 
buckling load with standard deviation of E n . The response shows a linear 
behavior. The plate with all edges CC is more affected by change in input 
randomness compared to a plate with all edges SS. The plate shows least 
sensitiveness to input dispersion for the c/d= 1.225. There is a marked 
increase in sensitiveness of the plate response for c/d= 0.4 compared to 
c/d= 1.225. The [0°/90 o /0°/9O o ] laminate is affected the least compared to 
angle-ply laminate for all edges SS due to change in SD of En. 
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4.7.2 Effect of transverse elastic modulus E u 

Fig 4.33 shows the variation of mean of the buckling load with standard 
deviation of E^. Fig 4.34 shows the variation of standard deviation of 
buckling load with standard deviation of E n . The general behavior is 
almost linear, but unlike the case E n , here the slope of both mean and 
standard deviation of normalized critical load curves are decreasing trend. 
The effect of random variable E n is less compared to random variable E n . 
As observed in the case of Eu here is again a marked higher sensitiveness 
for plates with c/d= 2.5 as compared to plate with c/d= 1.225. All edges 
CC plate is more affected by dispersion in E Vi than all edges SS. The 
[0 o /9O o /0°/9O°] laminate is affected lesser than [45°/ — 45°/45 0 / — 45°] laminate 
for all edges SS. 

4.7.3 Effect of major Poisson’s ratio v\ 2 

Fig 4.35 shows the variation of mean of the buckling load with standard 
deviation of i / 12 . Fig 4.36 shows the variation of standard deviation of 
buckling load with standard deviation of vu . The general nature of the 
curves are mostly linear. The plate with all edges CC is more affected than 
plate with all edges SS. From the comparison of all variables it has found 
that the variable v \2 is least effective compared to random variables Eu, 
E \2 and G\ 2 ’ 

4.7.4 Effect of rigidity modulus G\ 2 

Fig 4.37 shows the variation of mean of the buckling load with standard 
deviation of Gn> Fig 4.38 shows the variation of standard deviation of 
buckling load with standard deviation of (?i 2 * 

Here too the characteristics are mostly linear. The plate with all edges 
CC is more affected than plate with all edges SS. Next to the random 
variable E lu the G n is more effective for buckling characteristics. 



70 


4.7.5 Combined effect of longitudinal modulus E n and rigid- 
ity modulus G 12 

Fig 4.39 shows the variation of mean of the buckling load with standard 
deviation of E n and G 12 . Fig 4.40 shows the variat ion of standard deviation 
of buckling load with standard deviation of E n and G 12 

. The [0°/90°/0°/ 90°] laminate is least affected compared to [45°/ —45° /— 45° /45°] 
and [45°/— 45°/ 45°/— 45°] laminates with all edges SS. The plate with all 
edges CC is more affected than plate with all edges SS. The general be- 
havior is almost linear, unlike the cases of random variables En, F 12 , v n , 

Gn, here the slope of mean and standard deviation curves are increasing 
trend. 


4.8 Normalized Mean and SD of critical loads for 
various cut-out shapes for antisymmetric angle- 
ply laminate 

The results in Table 4.4 shows the variation of Mean of normalized critical 
loads with variation of SD. Table 4.5 shows the variation of SD of nor- 
malized critical loads with variation of SD. From the above tables, for a 
given area of the cut-outs, the buckling load for a laminate with elliptical 
cut-out is lower than circular and rectangular cut-outs. The buckling load 
with circular cut-out is highest. 
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Table 4.4: Comparison of normalized (Mean) critical loads for various cut-outs 


MATERIAL 

PROPERTY 

a/n 

RATIO 

E 

LLIPTICAL - CIT 

TOUT 

RECTANGULAR - CUTOUT 

CIRCULAR 

c/d=2.5 

c/d =1.225 

c/d=0.4 

c/d=1.15 

c/d=0.96 

c/d=0.87 

a/d=1.58 

En 

0 

1.4492395 

1.4240004 

1.4421920 

15621620 

1596388 

15682254 

1.6495797 

5 

1.5053510 

1.4926805 

1.4980849 

1.6224700 

1.6581648 

1.691938 

1.7116374 

10 

15609732 

1.5444479 

1.5534758 

1.6931614 

1.7193275 

1.6896032 

1.7732239 

15 

1.6194321 

1.5849438 

1.6086270 

1.7015999 

1.7801675 

1.7497385 

1.8346228 

20 

1.6711078 

1.6376066 

1.6631463 

1.8002852 

1.8402538 

1.8091720 

1.9860838 

E u 

0 

1.4492395 

1.4240004 

1.4421920 

1.5621620 

15963880 

1.5682254 

1.6495797 

5 

1.4513068 

1.4260352 

1.4443664 

15641722 

1.6184553 

1.5703854 

1.6532104 

10 

1.4533712 

1.4280678 

1.4465381 

15661797 

1.6355194 

15725422 

1.6568033 

15 

1.4554226 

1.4300868 

1.4486954 

15681739 

1.6500197 

15746845 

1.6603673 

20 

1.4574935 

1.4321247 

1.4508732 

15701863 

1.6716710 

1.57 68464 

1.6639507 

o 

1-1 

0 

1.4492395 

1.4240004 

1.442192 

15621620 

15963880 

15682254 

1.6495797 

5 

1.4569753 

1.4317882 

1.4488694 

15699273 

1.6041548 

1.5754868 

1.6603959 

10 

1.4626034 

1.4391237 

1.4551450 

15772181 

1.6114264 

15823038 

1.6660424 

15 

1.4687346 

1.4460861 

1.4610682 

15840928 

1.6182656 

15887318 

1.6736924 

20 

1.4745546 

1.4527068 

1.4666797 

15887056 

ISiB 

15948189 

1.6810122 

Dl2 

0 

1,4492395 

1.4240004 

1.4421920 


1 



5 

1.4493346 

1.4240810 

1.4422972 


1.5964704 

15684444 

1.6500954 

10 

1.4494208 

1.4241702 

1.4424129 

15632266 


156865174 

1.6506239 

15 

1.4495270 

1.4242630 

1.4425388 

15638256 

15966682 

15689761 

1.6511625 

20 

1.4492395 

1.4243495 

1.4426755 

15643357 

15967843 


1.6520514 

E u &G 12 

0 

1.4492395 

1.4240004 

1.4421920 

j 

15621620 

15983880 

E2H 

1.6495797 

5 

1.5145100 

1.4878887 

1.5070596 

1.6326538 

1.6571546 

1.6389737 

1.7446564 

10 

1374315 

1.5516605 

1.5718314 

1.7030315 

1.7403914 

1.7096150 

1.7950891 

15 

1.644764 


1.6365063 

i 

1.7732944 

1.8122135 

1.7801477 

1.8676808 

20 

■ 

1.7171183 

jgggjjjl 

1 1.7083344 

1.8518997 

1.8942260 

1.8584834 

1.9490373 
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Table 4.5: Comparison of normalized (SD) critical loads for various cut-outs 


MATERIAL 

a/p. 

E 

LLIPTICAL - CUTOU 

T 

RECTANGULAR - CUTOUT 

CIRCULAR 

PTJT . Of FT 

PROPERTY 

RATIO 

c/d=2.5 

c/d=1.225 

c/d=0.4 

c/d=1.15 

c/d=0.96 

c/d=0.87 

a/d =1.58 


0 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 


5 

0.035745 

0.0405153 

0.0356054 

0.0384147 

0.0407634 

0.0388382 

0.0395384 

E„ 

10 

0.071095 

0.0806950 

0.0708079 

0.0811873 

0,0810275 

0.0772276 

0.0786960 


15 

0.108770 

0.1192131 

0.1057207 

0.1139517 

0.1207210 

0.1152915 

0.1175915 


20 

0.140840 

0.1555576 

01402550 

0.1511025 

0.1557256 

0.1529321 

0.1678188 


0 

0.000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 


5 

0.001317 

0.0012967 

0.0013854 

0.0013987 

0.0013168 

0.0013979 

0.0023064 

En 

10 

0.002632 

0.0025914 

0.0027686 

0.0027990 

0.0026317 

0.0027930 

0.0045980 


15 

0.003935 

0.0038735 

0.0041389 

0.0042089 

0.0039393 

0.0041819 

0.0068564 


20 

0.005257 

0.0051745 

0.0055292 

0.0056708 

0.0054247 

0.0055987 

0.0091372 


0 

0.000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 


5 

0.004162 

0.0049255 

0.0042261 

0.0049143 

0.0050172 

0.0045956 

0.0068186 

G 12 

10 

0.008409 

0.0095218 

0.0081467 

0.0094676 

0.0098529 

0.0088526 

0.0103780 


15 

0.012197 

0.0138287 

0.0118016 

0.0137076 

0.014027 6 

0.0128168 

0.0151268 


20 

0.015753 

0.0178806 

0.0152248 

0.0163800 

0.0182987 

0.0165264 

0.0196296 


0 

0.000000 

f 

0.0000000 

0.0000000 

0.000000 

0.0000000 

0.0000000 

0.0000000 


5 

0.000055 

0.0000517 

0.0000677 

0.0000777 

0.0000741 

0.0000698 

0.0003293 


10 

0.000118 

0.0001111 

0.0001433 

0.0001467 

0.0001281 

0.0001503 

0.0006673 


15 

0.000189 

0.0001717 

0.0002275 

0.0002075 

0.0001975 

0.0002751 

0.0010158 


20 

0.000269 

0.0002383 

0.0003197 

0.0003481 

0.0002450 

0.0003427 

0.0013265 


0 

0.000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

0.0000000 


5 

| 0.040772 

0.0402224 

0.0404751 

0.0441455 

0.0477002 

0.0441413 

0.0619492 

E u &G 12 

10 

0.083770 

0.0805695 

0.0814567 

0.0883622 

0.0902348 

0.0883475 

0.0914177 


15 

0.122494 

0.1208631 

0.1216024 

0.1326824 

0.1354655 

0.1326209 

0.1371957 


20 

0.159243 

0.1564597 

0.1581719 

0.1721966 

0.175937 

0.1725041 

0.1779614 













































Chapter 5 


Conclusions 


In this thesis, a statistical study of the buckling behavior of composite 
laminated plates with different configurations has been carried out with the 
help of Monte Carlo simulation, finite element method and commercially 
available NASTRAN software. The following main conclusions have been 
drawn from this study. 

1. The longitudinal modulus En and the rigidity modulus G 12 are the 
most critical of the material properties. They considerably influence 
the buckling characteristics of the composite plate. 

2. As the AR increases the mean of the buckling load of the plate de- 
creases to the variation in the input material property characteristics. 

3. The normalized critical load characteristics of the plate is affected by 
the edge conditions. The ply-orientation also affects the character- 
istics but its effect is quite small compared to the effect of the edge 
conditions. 

4. The normalized critical load varies almost linearly with different ARs, 
boundary conditions and ply-orientations. 

5. The CC plate is more sensitive to the changes in material properties 
compared to SS plate. 
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6. For a give area of the cut-outs, the buckling load for a laminate with 
elliptical cut-out is lower than circular and rectangular cut-outs. The 
buckling load with circular cut-out is the highest. 


5.1 Scope for further work 

Having worked with this interesting problem for long the author would like 
to suggest the following for future work. 

1. The effect of dispersion in the random material properties on the 
stresses and dynamics of the plate can be studied. 

2. The random material properties were taken to have Normal distribu- 
tion over lots and was assumed not to vary over the domain of the 
plate. The problem may be pursued with the material properties hav- 
ing different probability distributions over the domain of the plate as 
well as varying in lots. 

3. Studies to be carried out with bi-axial loading, plate geometry being 
circular, elliptical etc. The boundary conditions may be fixed, SS, free 
and in combinations of different boundary conditions. The parameters 
thickness and ply-orientation, etc. may be taken as random in nature. 

4. Studies to be carried out with different cut-out shapes other than 
regular boundaries with random material properties. 
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